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The  characteristic  ‘N’  shaped  pressure  waves  recorded  at  ground  level  below  the  flight  path  of  supersonic 
aircraft  have  compressive  regions  very  much  less  steep  than  predicted  by  theory,  with  a consequent  reduction 
in  the  subjective  annoyance.  Attempts  have  been  made  to  explain  this  phenomenon  in  terms  of  the  effect  of 
atmospheric  turbulence.  This  Report  is  concerned  with  the  particular  proposal  that  since  on  passing  through 
turbulence  a shock  will  tend  to  become  thicker  for  part  of  the  time  and  to  steepen  for  the  remainder  then  if  the 
thickening  process  operates  faster  than  the  steepening  a net  increase  in  thickness  will  result. 

The  problem  is  simulated  in  one-dimension  by  a variation  of  the  ‘piston  problem’  in  which  pressure  waves 
representing  turbulence  are  made  to  pass  through  a shock  wave.  The  governing  partial  differential  equations 
are  solved  using  a finite  difference  technique. 

The  results  of  the  calculations  show  that  a shock  wave  tends  to  become  thicker  on  passing  through  a 
disturbance  wave  having  equal  positive  and  negative  parts.  This  is  consistent  with  the  proposed  explanation  of 
shock-thickening  although  the  thickening  calculated  is  small  compared  to  the  initial  shock  thickness,  for  the 
disturbance  waves  used.  Distortion  of  the  disturbance  waves  with  time  prevented  a wider  study  of  the  effects  of 
varying  amplitude  and  wavelength  being  made  so  that  precise  quantitative  conclusions  could  not  be  drawn  but 
it  is  suggested  that  the  process  investigated  is  unlikely  to  be  powerful  enough  to  produce  the  shock  thickening 
measured  experimentally  in  sonic  boom  research.  , 
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The  leading  and  trailing  shock  waves  generated  hy  an  aircraft  Hying  at  supersonic  speed  develop  with 
increasing  distance  from  the  aircraft  into  a characteristic  'N'  shaped  pressure  wave.  An  example  of  such  a 
wave,  generated  by  Concorde  and  reported  by  Holbeehe,  ct  at.,  is  given  in  Fig.  1 . This  example  is  typical  of 
such  waves  when  recorded  at  or  near  ground  level  below  the  flight  path  of  a high-flying  aircraft  in  that  while  the 
expansion  part  of  the  wave  conforms  with  theoretical  predictions,  the  compressive  regions  are  very  much 
thicker  than  expected.  The  factor  by  which  the  compressive  regions  are  thickened  can  be  as  high  as  10'.  The 
phenomenon  of  shock  thickening  is  of  considerable  practical  significance  in  that  it  is  found  to  reduce  the 
subjective  annoyance  resulting  from  the  sonic  boom. 

Atmospheric  turbulence  is  one  possible  cause  of  the  shock  thickening*  and  experiments  by  Bauer'  have 
demonstrated  some  thickening  with  the  introduction  of  artificially  produced  turbulence,  although  the  scale  of 
thickening  measured  in  the  experiments  was  at  least  an  order  of  magnitude  less  than  that  commonly  recorded 
under  the  flight  path  of  high-flying  supersonic  aircraft.  Several  investigators,  for  example,  Howe4  and  Plotkin 
and  George'  have  sought  an  explanation  of  the  thickening  in  the  scattering  of  waves  on  passing  through 
turbulence.  However,  in  a review  of  theories  based  on  wave-scattering  mechanisms,  Ffowcs  Williams  and 
Howeh  have  concluded  that  in  practice  a wave  may  be  thickened  in  this  way  by  a factor  of  2 at  the  most. 

An  alternative  explanation  of  shock  thickening  due  to  turbulence  has  been  proposed  by  Broadbent  . He 
argues  that  as  the  shock  enters  atmospheric  turbulence  the  equilibrium  which  may  have  existed  between  the 
nonlinear  steepening  and  the  thickening  due  to  diffusion  will  be  disrupted  by  large  local  disturbances  that  will 
be  repeated  in  a random  fashion  while  the  turbulence  persists.  Thus  each  part  of  the  shock  front  will  be  too 
thick  for  the  ambient  conditions  during  roughly  half  the  time,  when  the  nonlinear  effect  will  cause  it  to  steepen, 
and  for  the  remaining  time  be  too  thin,  when  diffusion  effects  will  tend  to  thicken  it.  The  two  opposing  effects, 
however,  operate  at  different  rates,  the  variation  due  to  nonlinear  steepening  depending  on  time  while  that 
due  to  diffusion  depends  on  the  square  root  of  time.  Thus,  in  small  intervals  of  time,  the  diffusion  process  will 
make  more  progress  than  the  steepening  effect  so  that  the  shock  will  become  progressively  thicker.  It  should  be 
possible  to  represent  the  essential  features  of  this  problem  by  considering  a variation  of  the  piston  problem  in 
which  a shock  wave  travels  down  a tube  and  passes  through  a pressure  wave  or  a series  of  waves  travelling  in 
the  opposite  direction.  This  Report  describes  a numerical  study  of  such  a representation  in  which  the 
governing  equations  are  solved,  using  a finite  difference  method,  and  the  thickening  of  the  shock  wave  by 
pressure  waves  of  differing  waveforms,  amplitudes  and  wavelengths  is  calculated.  Computational  considera- 
tions imposed  an  upper  limit  on  the  amplitude  and  wavelength  of  the  waves  representing  turbulence  and 
restricted  the  number  of  interactions  in  any  one  calculation  to  one  or  two.  However,  it  is  thought  that  the  scope 
of  the  calculations  is  sufficiently  wide  for  the  results  to  provide  some  insight  into  the  validity  of  the  suggested 
explanation  of  shock  thickening. 

The  simulation  of  the  problem  is  described  in  Section  2,  the  development  and  solution  of  the  finite  difference 
equations  is  explained  in  Section  3,  the  calculations  performed  are  described  in  Section  4 and  the  results  of  the 
calculations  are  given  in  Section  5 and  discussed  in  Section  6. 

The  results  of  the  calculations  show  that  some  shock  thickening  is  possible  but  the  computational  restraints 
prevent  an  estimation  of  the  possible  thickening  due  to  this  process  in  a real  case  of  a shock  wave  moving 
through  turbulence.  However,  it  is  inferred  that  on  passing  through  a series  of  similar  waves  representing 
turbulence  the  shock  wave  would  attain  a time-average  asymptotic  thickness  and  that  the  total  thickening  is 
unlikely  to  be  many  times  greater  than  that  produced  by  the  passage  of  a single  wave.  It  is  therefore  concluded 
that  the  process  considered  is  unlikely  to  be  the  principal  cause  of  the  observed  shock  thickening. 


2.  Simulation  of  the  Problem 

If  the  phenonemon  of  shock  thickening  on  passing  through  turbulence  is  indeed  due  to  the  different  rates  at 
which  steepening  and  thickening  operate  then  the  essential  features  of  the  problem  should  be  present  in  one 
dimersion  in  the  variation  of  the  ‘piston  problem-  illustrated  in  Fig.  2.  At  one  end  of  a long  tube  containing 
fluid  there  is  a piston  moving  with  a uniform  velocity  and  sending  a shock  wave  down  the  tube  while  at  the 
other  end  a second  piston  is  moving  in  an  oscillatory  manner  and  is  sending  a series  of  compression  and 
expansion  waves  in  the  opposite  direction  to  meet  and  pass  through  the  oncoming  shock  wave.  If  the  motion  of 
the  latter  piston  is  such  that  the  waves  it  produces  have  relatively  large  amplitudes  and  small  gradients, 
compared  to  the  shock  wave,  then  the  waves  will  resemble  turbulence  in  their  tendency  to  swamp  the  shock  by 


* Real  gas  effects'  have  also  been  proposed. 
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an  imposed  velocity  distribution  that  does  not  contain  any  sharp  fronts  such  as  the  shock  itself.  Real  turbulence 
will  have  other  three-dimensional  characteristics  that  may  be  important  in  the  interaction  with  the  shock  (e.g. 
vortieity),  but  the  simple  waves  used  in  the  present  paper  are  thought  to  be  sufficient  to  test  the  theory  outlined 
above,  which  depends  only  on  disturbing  the  equilibrium  of  the  shock. 

Throughout  this  Report  the  shock  will  be  regarded  as  having  originated  from  the  left  and  the  waves 
representing  turbulence  from  the  right,  as  in  Fig.  2. 

The  problem  may  be  simplified  if  the  pistons  are  regarded  as  being  beyond  the  length  of  tube  along  which  the 
calculation  is  to  be  carried  out.  The  effect  of  the  piston  at  the  left  may  then  be  represented  by  a constant 
particle  velocity  at  the  left-hand  boundary  and  the  effect  of  the  oscillating  piston  by  a particle  velocity  varying 
in  a suitably  prescribed  manner  at  the  right-hand  boundary.  The  boundaries  referred  to  are  the  two  ends  of  the 
length  of  tube  over  which  the  calculation  is  to  be  carried  out.  In  this  way  difficulties  in  the  numerical  solution 
resulting  from  having  a varying  length  of  tube  to  consider  are  avoided. 

Fig.  3 illustrates  an  example  of  the  variation  of  particle  velocity,  denoted  by  v,  with  distance  measured  along 
the  axis  of  the  tube,  denoted  by  x , before  the  two  waves  moving  in  opposite  directions  have  met.  The  wave 
representing  turbulence  could  be  repeated  to  the  right  indefinitely  although,  as  already  noted,  the  computa- 
tions reported  here  were  limited  to  one  or  two  cycles  only. 

A numerical  method  is  required  that  will  enable  the  wave  representing  a shock  to  be  developed  to  its 
asymptotic  form  and  allow  the  wave  or  waves  representing  turbulence  to  move  to  the  left  with  little  or  no 
distortion.  If  the  local  particle  velocities  can  then  be  calculated  as  the  shock  passes  through  and  beyond  the 
simulated  turbulence  it  will  be  possible  to  calculate  the  effect  of  the  turbulence  upon  the  shock  thickness. 

3 Numerical  Method 

3.1.  Basic  Equations 

l.ighthilIKhas  given  the  equations  of  plane  sound  waves  of  finite  amplitude  with  linearized  diffusion  terms  as 


dv  dv  2 da  „ d2V 

+ 1>  + a = S -> 

dt  dx  y—\dx  dx~ 

(1) 

da  da  y - 1 dv  . 

dt  dx  2 dx 

(2) 

where  t is  time,  a is  the  local  speed  of  sound,  y is  the  ratio  of  specific  heats  and  S is  the  ‘diffusivity  of  sound', 
that  is  the  combination  of  different  diffusivities  which  affect  attenuation  of  sound  waves.  The  left-hand  sides  of 
equations  (1)  and  (2)  are  the  exact  forms  of  the  equations  of  sound  waves  of  finite  amplitude  under 
thermodynamically  reversible  conditions  and  the  right-hand  sides  are  a ‘linearized'  approximation  to  the 
effects  of  diffusion. 

If  the  variables  in  equations  (1)  and  (2)  are  made  non-dimensional  with  respect  to  the  speed  of  sound  in 
undisturbed  fluid,  denoted  by  a„,  and  to  some  time  yet  to  be  defined,  denoted  by  so  that 


(3) 


(41 

Now  t„  may  be  chosen  to  make  D = 1,  that  is  /„  = 8/al,  so  that  equation  (3)  becomes 


v = a0V,  a=a0(l+A),  t = t0T, 
x = a0toX  and  S = alt»D 


then  equations  (1)  and  (2)  may  be  written,  respectively,  as 


and 


dV  ,,dV  2 , , M 

Jf+VM+^i(1+A)M=DvF 


dA  ,dA  y-1  , „ dV  n 

tf+VM  + — (l+A)M  = °- 


For  a numerical  solution  it  is  convenient  to  transform  the  x co-ordinate  further  so  that  the  shock  remains 
nearly  stationary  in  the  new  co-ordinate.  Writing  £ = X-riT  and  r=T  equations  (4)  and  (5)  become, 
respectively. 


and 
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(6) 


(7) 


To  hold  the  shock  nearly  stationary,  in  terms  of  (,  17  must  have  a value  of  approximately  I but  this  value  may  be 
modified  with  experience. 


3.2.  Finite  Difference  Equations 

To  establish  a finite  difference  scheme  the  variables  r and  £ are  treated  as  rectangular  co-ordinates  and  a 
mesh  with  equal  elements  At  and  A£  is  adopted.  Suffices  m and  n are  used  to  denote  position,  thus  Vm  n is  the 
value  of  V at  r = m At  and  £ = n A£.  A difference  scheme  that  is  ‘central’  with  respect  to  mesh  intervals  in  both 
t and  {;  is  adopted. 

Considering  equation  (6)  first,  we  write: 
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(8) 


(9) 


(10) 


(11) 


(12) 


Substituting  the  above  in  equation  (6),  which  is  equivalent  to  evaluating  equation  (6)  at  t = {m 
£ = (n  ~2)A£,  gives: 


where 


and 


P nAm-n  QnAmM  - 1 rn 

(13) 

Pn  = Wn  + V„  4-  zn, 

(14) 

qn  = H’„  - yn  + zn 

(15) 

\ )Am  1 <M  4*  ( wn  4-  yn  zn  )Am  - j n 1 

(16) 

with 


a step-by-step  solution  will  involve  an  iteration  process.  The  coefficients  of  equations  (13)  and  (26)  are  first 
evaluated  after  estimating  the  variables  by  some  means  and  then  the  equations  are  solved  in  their  linear  form 
to  obtain  new  values  of  the  variables.  These  are  then  used  to  improve  the  values  of  the  coefficients  and  the  cycle- 
repeated  until  convergence  is  achieved.  The  linear  form  adopted  in  equations  (13)  and  ( 16)  is  to  some  extent  a 
matter  of  choice,  in  that  with  a nonlinear  product  of  two  variables  either  variable  may  be  retained  explicitly 
while  the  other  is  relegated  to  the  coefficients,  but  experience  shows  that  stability  is  improved  by  keeping  the 
derivatives  explicit. 

3.3.  Solution  of  Finite  Difference  Equations 

In  order  to  solve  equations  (13)  and  (26)  certain  boundary  conditions  must  be  specified;  these  are  given  by 
the  prescribed  particle  velocities  at  each  end  of  the  length  of  tube  through  which  the  equations  are  to  be  solved 
and,  initially,  by  an  assumed  particle-velocity  distribution  along  the  tube  at  r = 0.  That  is,  if  the  equations  are 
to  be  solved  from  f = 0 to  £ = ,\'A£  then  Vm  V',„N  and  A,„v  must  be  specified  for  all  m and  V()  n and  A„n  for 
all  n.  A„,..\  is  found  from  Vm  N using  the  equation  for  local  acoustic  speed  (which  holds  in  the  absence  of 
diffusion) 

a = a0  + s(y- l)u,  (31) 

which  may  be  written  as 

Am.n=ky-l)VmM.  (32) 

Given  VUn  it  is  assumed  that  VY„  = V'„„  for  all  n.  Ax„  is  then  calculated  from  equation  (13),  working  down 
from  A i .v  and  re-writing  equation  (13)  as 

Am,„ - 1 = (r„  - p„A»,.n )/<?-•  (33) 

Now  equation  (26)  is  solved  for  VX  n by  the  resolution  of  a tri-diagonal  matrix,  using  the  values  of  Ai,„  just 
calculated  and  the  assumed  values  of  V',  „ to  evaluate  the  coefficients.  The  new  values  of  V',  „ are  then  used 
with  equation  (33)  to  find  improved  values  of  A, and  thus  an  iterative  procedure  is  established.  When 
successive  approximations  to  Vin  differ  by  less  than  a prescribed  tolerance  the  calculation  is  advanced  to 
compute  14, „ and  A2  n and  so  on  to  higher  values  of  m.  After  the  first  step  in  m is  completed,  linear 
extrapolation  is  used  for  the  first  estimate  of  V,„  „ at  the  end  of  each  subsequent  step.  The  first  and  last  rows  of 
the  tri-diagonal  matrix  are  written  as 

A,  VmA  +C|  Vm.j  = dx -a,  ,,  (34) 

and 

a.v  I Vm.s  : , VVv  i = tf.v  i - Cv  I V'm.s,  (35) 

respectively,  thus  introducing  the  boundary  conditions  and  V'm..v. 

The  above  procedure  has  been  programmed  for  a digital  computer  in  FORTRAN  language;  a standard  I CL 
subroutine  was  used  for  the  solution  of  the  tri-diagonal  matrix. 

4.  Calculations  Performed 

4.1.  The  Asymptotic  Shock 


In  order  to  study  the  effect  of  the  simulated  turbulence  upon  the  shock  it  was  first  necessary  to  establish  the 
asymptotic  shock  in  the  absence  of  disturbances  and  to  determine  its  thickness.  Lighthill*.  has  given  the 
thickness  of  an  asymptotic  shock  as 


where  h is  defined  as  the  distance  in  which  v changes  from  0-05uo  to  0-95t>0  and  v,>  is  the  particle  velocity  on 
the  compressive  side  of  the  shock.  The  non-dimensional  shock  thickness  is  therefore  given  by 
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. t 
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H = 


12 

(y  + 1 ) k'n 


where  V„  is  the  non-dimensional  particle  velocity  corresponding  to  v„.  A value  of  ()•  1 was  assigned  to  V7,,;  this 
is  large  for  a sonic  boom  but  one  of  the  problems  of  the  numerical  calculation  was  to  obtain  significant  effects 
from  a computer  program  without  making  extravagant  demands  on  working  storage  space.  Equation  (37) 
shows  that  a shock  thickness  of  approximately  50  units  in  f may  be  expected  for  a value  of  V(,  of  0- 1 . To  start 
the  calculation  a linear  approximation  to  the  shock  was  used,  with  a thickness  of  40  units  in  £.  Thus  V„  „ was 
equal  to  0- 1 from  £ = 0 to  some  higher  value  of  £,  then  decreased  linearly  with  £ to  zero  over  40  units  in  £ and 
remained  at  zero  to  £ = NA£  A„.n  was  estimated  from  V„  n for  all  n,  using  equation  (32).  The  calculation  was 
then  advanced  in  time  as  described  above,  keeping  V,,,  „ = 0- 1 and  Vm  N = 0,  until  the  slope  and  thickness  of  the 
shock  remained  constant  with  time. 

In  order  to  check  that  the  calculated  asymptotic  shock  form  is  the  same  irrespective  of  whether  the  initial 
approximation  was  ‘too-steep’  or  'too-thick’,  a second  computation  was  made  similar  to  the  above  but  with  the 
initial  linear  approximation  extending  over  60  units  in  £.  This  did  in  fact  confirm  that  the  asymptotic  shape  was 
independent  of  the  starting  assumptions. 

Two  further  computations  were  made,  starting  from  shock  forms  that  were  obtained  by  scaling  the 
asymptotic  shock  with  respect  to  f ; one  calculation  was  started  from  a ‘too-thick’  and  the  other  from  a 
‘too-steep’  shock.  The  asymptotic  shock  form  was  then  developed  in  each  case,  thus  allowing  a comparison  to 
be  made  between  the  rates  at  which  the  shock  thickened  and  steepened. 

4.2.  Waves  Representing  Turbulence 

Having  established  the  asymptotic  shock  form,  waves  representing  turbulence  were  introduced  by  varying, 
with  increeasing  r,  the  boundary  condition  at  £ = 7VA£,  that  is  at  the  right-hand  end  of  the  region  of 
computation.  Calculations  were  carried  out  with  waves  of  three  different  forms  to  be  referred  to  as  types  (A), 
(B)  and  (C),  the  respective  particle  velocities  at  the  boundary  being  given  by 

Vm.N  =s{sin  [2(1  + i7)ttt/A]},  (38) 

Vm.N  = s{l  - COS  [2(  1 + 17)777-/ A ]}/2  (39) 

and 

Vm.n  - *{1  - cos  [2(1  + t i)ttt/\]}  sin  [2(1  + rj)T7r/A]/{[l  - cos  (27r/3)]sin  (2ir/3)}.  (40) 

Here  s defines  the  amplitude  (although  its  sign  is  significant)  and  A the  wavelength  of  the  waves  as  they  were 
introduced  at  the  boundary.  Values  of  A„ -N  corresponding  to  Vm,N  were  computed  from  equation  (32).  The 
variation  of  the  right-hand  boundary  condition,  as  given  by  equations  (38),  (39)  and  (40)  is  illustrated  in  Fig.  4. 
The  waves  produced  travel  to  the  left  with  a velocity,  in  the  transformed  co-ordinate  £ of  approximately  2 
units  per  unit  advance  in  r. 

Values  of  s from  -1-0  to  1-0  and  of  A from  250  to  1000  were  used;  details  are  to  be  found  in  Section  5. 
Practical  considerations  imposed  an  upper  limit  on  A such  that  the  values  used,  though  large  compared  to  the 
shock  thickness,  were  small  compared  to  the  probable  wavelengths  of  turbulence  in  the  atmosphere. 

Two  additional  calculations  were  made  with  wave  forms  of  type  B in  which  at  the  start  of  the  computation 
the  asymptotic  shock  form  was  replaced  by  one  about  10  per  cent  thicker.  The  object  of  these  two  calculations 
was  to  investigate  the  effect  of  waves  representing  turbulence  upon  a shock  initially  having  a strong  tendency 
to  revert  to  its  asymptotic  thickness. 

4.3.  Boundary  Condition  at  f = 0 

As  the  wave  representing  turbulence  must  pass  through  and  clear  of  the  shock  before  its  effect  on  shock 
thickness  can  be  calculated,  a very  long  region  of  calculation  is  necessary  to  accommodate  the  wave  unless 
some  means  can  be  found  of  ‘absorbing’  it,  without  ‘reflections’  occurring,  at  the  left-hand  boundary  of  the 
calculation.  Computational  considerations  impose  a practical  limit  on  the  length  of  the  region  considered  and 
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the  difficulty  is  aggravated  by  the  tendency  of  the  waves  to  increase  in  wavelength  with  time  and  to  be  followed 
by  'tails’  of  waves  of  decreasing  amplitude,  both  these  effects  being  inherent  in  the  numerical  technique 
employed.  It  was  therefore  necessary  to  devise  a means  of  modifying  the  boundary  condition  at  £ = 0 from  the 
initial  condition  Vm  „ = Vit  in  order  to  absorb  the  turbulent  wave.  As  the  wave  representing  turbulence  travels 
to  the  left  with  a velocity  of  approximately  2 units  in  ^ for  each  unit  increase  in  r it  is  possible  to  estimate  a 
value  of  V at  a given  point  from  a value  or  values  to  the  right  at  previous  steps  in  r.  Thus,  for  example,  a rough 
approximation  is  given  for  Vm  „ by 

Vm,n**Vm  (4H 

assuming  that  the  step  lengths  At  and  Af  are  equal.  By  this  means  it  was  possible  to  estimate  values  for  the 
left-hand  boundary  conditions  approximating  to  the  values  of  V and  A that  would  have  occurred  at  the  same 
point  but  with  the  region  of  calculation  extending  infinitely  far  to  the  left.  In  order  to  make  a correction  for  the 
effects  of  wave  distortion  and  variation  in  wave  speed,  second  order  approximations  for  the  boundary 
conditions  were  adopted,  thus 

Kb.o  = 2 Vm_i,2— V^-i.4,  (42) 

where  again  it  is  assumed  that  Ar  = A£  Where  Ar  and  A£  were  not  equal,  equation  (42)  was  modified  as 
appropriate.  For  example  with  At  = 10  and  A£  = 4 

Vm, 0=  2 Vm_!  5 — V,„_2,10  (43) 

was  used. 

As  a check  on  the  effectiveness  of  the  above  technique,  calculations  of  the  effect  of  turbulence  upon  shock 
thickness  were  made  with  the  shock  at  two  different  initial  positions.  The  results  of  the  two  calculations  were 
substantially  the  same,  demonstrating  that  the  effect  of  any  reflection  from  the  left-hand  boundary  condition 
was  negligible,  since  it  would  not  be  equal  in  the  two  calculations. 

4.4.  Calculations  of  Shock  Thickness 

In  order  to  calculate  the  thickness  of  the  shock,  H,  as  defined  in  Section  4. 1 , the  values  of  £ corresponding  to 
0 05  V()  and  0-95  V„  were  found  by  inverse  interpolation  from  the  values  of  V calculated  at  incremental  points 
in  t 

To  calculate  the  increase  in  H due  to  the  passage  of  the  wave  or  waves  representing  turbulence,  the  value  of 
H is  required  at  a value  of  t corresponding  to  the  instant  at  which  the  last  wave  has  just  cleared  the  shock;  any 
further  advance  in  t allows  the  shock  to  begin  the  process  of  reverting  to  its  asymptotic  thickness  and  the  effect 
on  H of  the  passage  of  the  waves  is  apparently  reduced.  However,  it  was  found  that  values  of  H computed  for 
successive  increments  in  r just  after  the  last  wave  representing  turbulence  had  cleared  the  shock,  fluctuated  to 
such  an  extent  that  it  was  impossible  to  estimate  H accurately  directly;  these  fluctuations  were  traced  to  the 
presence  of  a ‘tail'  of  minor  waves,  arising  from  the  numerical  method,  following  behind  the  wave  or  waves 
representing  turbulence.  The  calculation  was  therefore  advanced  in  t until  the  fluctuations  were  damped  out 
sufficiently  for  H to  be  estimated  with  confidence.  The  resulting  value  of  H was  then  corrected  to  allow  for  the 
amount  the  shock  would  have  changed  in  thickness  in  the  time  that  had  elapsed  since  the  wave  or  waves 
representing  turbulence  had  just  cleared  it,  using  the  values  of  H obtained  in  the  calculation  of  the  asymptotic 
shock  starting  from  too-steep’  or  'too-thick'  shapes. 

4.5.  Computation  Parameters 

After  numerical  experiments  to  establish  a good  compromise  between  computing  time  and  accuracy  the 
following  values  for  certain  parameters  were  used  in  the  majority  of  the  calculations:  N = 400,  At  = 10  and 
A£  = 4.  A value  of  1 -06  for  the  transformation  constant  r\  in  equations  (6)  and  (7)  was  found  to  hold  the  shock 
nearly  stationary  with  respect  to£  while  t was  advanced,  in  the  absence  of  simulated  turbulence.  Waves  of  type 
B had  the  effect  of  moving  the  shock  to  the  left  if  s was  negative  and  to  the  right  if  s was  positive.  Values  of  rj 
ranging  from  0-8  to  1 -3.  as  appropriate  to  a given  case,  were  therefore  used  to  hold  the  shock  near  the  centre, 
with  respect  to  £,  of  the  range  of  computation. 
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5.  Results  of  Calculations 


5.1.  The  Asymptotic  Shock 

The  calculated  asymptotic  shock  is  illustrated  in  Fig.  5.  The  thickness  of  the  shock,  as  defined  in  Section  4- 1 , 
is  48-6  units  in  f.  A similar  shock  shape  was  obtained  with  calculations  started  from  both  too-thin'  and 
•too-thick’  linear  approximations. 

The  result  of  calculations  made  starting  from  initial  shock-shapes  obtained  by  scaling  the  asymptotic  shock 
with  respect  to  f is  shown  in  Fig.  6;  SH  is  the  change  in  shock  thickness  corresponding  to  an  increment  St  in  r. 
It  may  be  seen  that  there  is  a slight  tendency  for  a too-thin  shock  to  thicken  at  a faster  rate  than  that  at  which  a 
too-thick  shock  steepens. 

5.2.  Effect  of  Waves  Representing  Turbulence 

5.2.  1.  Waves  of  Type  A.  Calculations  of  the  thickness,  H,  of  the  shock  after  one  wave  of  type  A had  passed 
through  it  showed  that  the  value  of  H was  fluctuating  with  successive  increments  in  r.  This  effect  was  traced  to 
the  presence  of  a succession  of  minor  waves  of  progressively  smaller  amplitudes,  following  behind  the  wave 
representing  turbulence.  These  minor  waves  were  found  to  be  present  to  some  extent  in  the  results  of  the 
calculations  for  all  three  types  of  wave  considered,  and  a method  of  overcoming  the  resultant  difficulty  in 
estimating  the  effect  on  shock  thickness  of  the  passage  of  the  wave  representing  turbulence  has  been  described 
in  Section  4.  However,  for  waves  of  type  A the  fluctuations  in  the  values  of  H were  of  such  magnitude  that  it 
was  impossible,  by  advancing  r,  to  obtain  an  accurate  value  for  H before  the  shock  had  reverted  to  its 
asymptotic  shape;  thus  no  useful  results  could  he  obtained  using  waveforms  of  type  A.  The  severity  of  the 
fluctuations  in  values  of  H for  successive  values  of  r was  probably  due  to  the  discontinuity  in  the  variation  of 
the  particle  velocity  at  the  abrupt  termination  of  the  initial  wave. 

5.2.2.  Waves  of  Type  B.  The  results  of  calculations  carried  out  with  turbulence  represented  by  waveforms  of 
type  B are  given  in  Table  1 ; AH  is  the  increment  in  shock  thickness  caused  by  the  passage  of  the  turbulence.  It 
can  be  seen  that  the  changes  in  H are  small  relative  to  the  asymptotic  value  for  H of  48-6  and  that,  in  general 

ti)  a negative  wave  increases  shock  thickness 

(ii)  a positive  wave  reduces  shock  thickness 

(iii)  the  influence  of  a wave  increases  with  amplitude  and  wavelength  and  has  a cumulative  effect  over  a 
number  of  cycles. 

However,  it  is  difficult  to  draw  any  quantitative  conclusions  because  of  the  tack  of  system  in  the  variation  of  AH 
with  the  waveform  parameters. 

In  order  to  find  an  explanation  for  the  inconsistent  behaviour  of  the  computed  values  of  AH  the  propagation 
of  the  waves  representing  turbulence  across  the  region  of  computation  has  been  examined.  Calculations  of  the 
development  of  these  waves  have  been  made  in  the  absence  of  a shock;  examples  of  the  results  of  these 
calculations  are  given  in  Figs.  7 and  8 where  the  variation  of  particle  velocity  with  r at  several  values  of  £ for  a 
wavelength  of  500  is  illustrated.  The  results  of  the  calculations  show  that  the  waves  undergo  distortion,  which 
becomes  more  severe  with  increasing  amplitude  and  decreasing  wavelength,  and  that  the  waves  for  positive 
values  of  s assume  different  shapes  from  their  counterparts  for  negative  values.  Furthermore,  because  of  the 
steepening  of  the  region  of  decreasing  velocity  in  some  cases,  for  example  s = TO  in  Fig.  7,  the  'tail'  of  minor 
waves  has  developed  to  such  an  extent  that  some  effect  on  shock  thickness  seems  probable.  However,  in  the 
majority  of  cases  the  distortion  has  not  produced  substantial  changes  in  either  the  amplitudes  or  the  lengths  of 
the  waves  and  it  seems  probable  that  if  the  analysis  is  confined  to  these  cases  then  the  results  should  have  some 
significance.  On  this  basis,  the  following  results  in  Table  1 have  been  bracketed  and  should  be  disregarded: 


A = 250,  s = 0-5  (1  and  2 cycles) 

A =500,  5=1-0. 

The  remaining  values  of  AH  have  been  plotted  against  s in  Fig.  9;  the  scatter  in  the  results  is  thought  to  be  due 
to  the  presence  of  the  minor  waves,  as  described  in  Section  4.4.  the  effect  of  which  may  not  have  been 
completely  eliminated.  Despite  the  scatter  it  is  apparent  that  there  is  a marked  tendency  for  a negative  wave  to 
have  a greater  effect  on  shock  thickness  than  a positive  wave  of  the  same  amplitude. 
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The  results  of  the  two  additional  calculations  that  were  made  with  the  initial  shock  having  a thickness  about 
10  per  cent  greater  than  that  of  the  asymptotic  shock  are  given  below: 


A 

S 

500 

-0-5 

-1-2 

500 

0-5 

-4-4 

— 

— 

-3-2 

The  third  of  the  above  three  cases  gives  the  steepening  of  the  initial  shock  during  the  same  advance  in  r as 
occurred  with  the  other  two  cases  but  in  the  absence  of  waves  representing  turbulence.  It  may  be  seen  that  in 
the  first  case,  although  the  shock  has  steepened  it  is,  after  the  passage  of  the  wave,  2-0  units  in  £ thicker  than 
the  third  case  at  the  corresponding  value  of  r.  The  value  of  2-0  may  be  compared  with  a similar  value  in  the 
corresponding  case  in  Table  1.  The  difference  between  the  second  and  third  cases  yields  a value  for  \H  of 
-T2,  which  may  be  compared  to  one  of  -T6  for  the  corresponding  case  in  Table  1.  Thus  it  seems  that  the 
thickening  due  to  the  tendency  of  the  shock  to  revert  to  its  asymptotic  thickness  may  be  added  to  the 
thickening  or  steepening  due  to  the  effect  of  the  wave  or  waves  representing  turbulence  to  give  the  net 
thickening.  If  this  is  in  fact  the  case  then  the  calculated  effects  of  the  waves  representing  turbulence  upon  the 
asymptotic  shock  must  also  include  the  opposing  effect  of  either  diffusion  or  non-linear  steepening,  although 
this  should  be  small  in  most  cases  since  the  rate  of  change  of  thickness  with  time  is,  of  course,  small  for 
thicknesses  close  to  the  asymptotic,  as  shown  in  Fig.  6. 

5.2.3.  Waves  of  Type  C.  As  with  waves  of  type  B,  it  was  found  that  the  waves  of  type  C t i ded  to  become 
distorted  as  they  moved  to  the  left  and  that  in  some  cases  the  wave  shape  at  the  position  of  tne  shock  was  so 
unlike  the  initial  shape  that  the  results  could  not  be  used  for  comparative  purposes.  However,  in  some  cases 
where  the  velocity  gradients  in  the  waves  were  sufficiently  small  the  waves  retained  their  essential  characteris- 
tics. Examples  of  such  waves  are  shown  in  Fig.  10  which  illustrates  the  variation  of  particle  velocity  with  r for 
A = 1000,  s=0-5  and  A = 1000,  s = -0-5  at  a station  corresponding  approximately  to  the  position  of  the 
i centre  of  the  shock  when  present,  although  Fig.  10  is  derived  from  calculations  made  in  the  absence  of  a shock. 

. i To  assist  in  comparing  the  shapes  of  the  two  waves  the  initial  waves,  that  is  the  waves  introduced  at  the  right 

hand  boundary,  have  been  added  to  Fig.  10.  It  may  be  seen  that  the  distortion  is  not  severe  and  that,  in 
particular,  the  positive  and  negative  parts  of  the  waves  have  remained  similar. 

Results  of  the  calculations  of  the  effect  on  shock  thickness  of  waves  of  type  C are  given  in  Table  2;  only  cases 
where  the  waves  representing  turbulence  were  not  severely  distorted  are  retained.  It  is  apparent  that  except 
for  one  case,  which  is  of  small  amplitude,  the  effects  of  the  negative  and  positive  parts  of  the  waves  do  not 
cancel  out  and  that  the  net  result  is  an  increase  in  shock  thickness;  this  is  consistent  with  the  results  of  Table  1. 
The  difference  between  the  values  of  H for  positive  and  negative  values  of  5 for  a wavelength  of  1000  can  be 
explained  by  the  effect  deduced  above  from  the  results  obtained  from  the  interaction  of  waves  of  type  B with  a 
shock  initially  too-thick.  If  the  positive  part  of  the  wave  passes  through  the  shock  first  its  steepening  effect  will 
be  opposed  throughout  by  the  tendency  of  the  shock  to  revert  to  its  asymptotic  shape  but  none  the  less  a net 
steepening  will  be  present  when  the  negative  part  of  the  wave  begins  its  passage  through  the  shock.  The 
thickening  effect  of  the  negative  part  of  the  wave  will  therefore  be  augmented  by  the  tendency  of  the  shock  to 
revert  to  its  asymptotic  thickness  until  such  time  as  the  asymptotic  thickness  is  achieved.  Thereafter  the  two 
effects  will  be  in  opposition  but  it  is  evident  that,  no  matter  whether  the  positive  part  of  the  wave  precedes  or 
follows  the  negative  part,  the  part  that  passes  through  the  wave  last  will  have  its  effect  on  the  shock  thickness 
increased.  The  results  in  Table  2 for  the  smaller  waves  conflict  with  this  hypothesis,  but  the  calculated  values  of 
H are  very  small  and  may  be  dominated  by  errors. 

6.  Discussion  of  Results 

In  the  preceding  Section  numerical  results  are  given  that  show  that  the  thickening  of  the  shock  caused  by  an 
interaction  with  a negative  wave  is  greater  than  the  steepening  produced  by  an  interaction  with  a positive  wave 
having  the  same  amplitude  and  wavelength.  This  effect  is  confirmed  by  the  results  for  a single  wave  having 
equal  positive  and  negative  parts,  and  it  is  consistent  with  the  proposed  explanation  of  shock  thickening  in 
which  it  is  suggested  that  the  two  opposing  effects,  diffusion  and  steepening,  operate  at  different  rates.  Some 
additional  support  to  the  suggestion  is  given  by  Fig.  6 which  shows  that  a shock  thinner  than  the  asymptotic 
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shock  thickens  at  a faster  rate  than  that  at  which  a too-thick  shock  steepens.  Other  results  are  given  that 
demonstrate  that  the  influence  of  a wave  increases  with  amplitude  and  wavelength  and  that  there  is  a 
cumulative  effect  on  thickening  over  a number  of  cycles,  although  calculations  were  made  for  one  and  two 
cycles  only. 

The  results  for  calculations  made  to  find  the  effect  of  waves  passing  through  a shock  initially  too-thick 
showed  that  the  steepening  process  goes  on  even  during  the  passage  through  the  shock  of  a wave  that  is  tending 
to  thicken  it,  the  net  effect  being  approximately  equal  to  the  sum  of  the  two  thickness  changes  calculated 
separately.  If  a shock  were  to  pass  through  a series  of  similar  waves  representing  turbulence  then  it  seems  likely 
that  a time-average  asymptotic  thickness  would  be  reached  when  the  rate  of  thickening  due  to  the  “turbulence- 
would  balance  that  due  to  the  tendency  of  the  shock  to  revert  to  its  undisturbed  asymptotic  thickness.  For 
example,  a succession  of  waves  with  A = 500  and  s = -0-5,  which  individually  would  thicken  the  initial  shock 
by  2-0  units  in  £ in  a time  of  approximately  250  units  in  r,  would  produce  an  asymptotic  thickness,  deduced 
from  Fig.  6,  of  about  52-5  units  in  £.  Thus,  in  this  case,  an  infinite  number  of  waves  would  produce  a thickening 
which  is  less  than  10  per  cent  of  the  initial  shock  thickness  and  only  twice  that  produced  by  a single  wave.  An 
increase  in  wavelength  of  the  turbulence  would  increase  the  effect  on  thickness  of  the  first  wave  but  the  effect 
of  subsequent  waves  would  soon  be  cancelled  out  by  the  rapidly  increasing  and  continuous  opposition  of  the 
steepening  effects,  which  would  have  longer  to  operate  during  each  interaction.  An  increase  in  amplitude 
would  also  increase  the  thickening  due  to  the  first  wave  but  the  second  wave  would  have  a greatly  reduced 
thickening  effect  since  the  steepening  tendency  would  be  much  increased  at  the  new  shock  thickness,  and 
equilibrium  would  be  reached  after  the  passage  of  fewer  waves  than  for  waves  of  smaller  amplitude.  It  seems 
unlikely,  therefore,  that  the  mechanism  investigated  here  could  be  the  cause  of  the  observed  shock  thickening. 

7.  Conclusions 

The  effect  of  turbulence  upon  the  thickness  of  a shock  wave  has  been  investigated  using  a numerical  method 
in  which  the  turbulence  was  represented  by  pressure  waves  which,  because  of  computational  restraints,  were 
not  more  than  two  in  number.  Distortion  of  the  waves  occurred  as  the  calculations  were  advanced  in  time  and 
resulted  in  waves  of  the  same  family  becoming  non-similar;  consequently  some  of  the  results  had  to  be 
discarded.  Disturbance  waves  having  initially  small  velocity  gradients  were  less  subject  to  distortion  and  the 
results  from  calculations  carried  out  with  such  waves  indicated  a tendency  for  the  shock  to  thicken  after  passing 
through  waves  having  equal  positive  and  negative  parts.  The  shock  thickening,  for  the  limited  range  of 
amplitudes  and  wavelengths  used  in  the  calculations,  was  small  compared  to  the  initial  shock  thickness  but  the 
trend  of  the  results  is  in  accordance  with  the  suggestion  that  on  passing  through  turbulence  a shock  wave  will 
tend  to  become  thicker  for  part  of  the  time  and  to  steepen  for  the  remainder,  with  the  thickening  process 
operating  at  the  faster  rate  so  that  a net  increase  in  thickness  results,  as  shown  by  Fig.  9.  Some  additional 
support  to  the  suggestion  is  given  by  Fig.  6 which  shows  that  a too-thin  shock  thickens  at  a faster  rate  than  that 
at  which  a too-thick  shock  steepens. 

Further  calculations  of  the  effect  of  a single  wave  upon  an  initially  too-thick  shock  showed  that  the  tendency 
for  the  shock  to  steepen  could  be  dominant  over  any  thickening  effect  due  to  the  disturbing  waves.  This 
suggested  that  on  passing  through  a series  of  similar  waves  representing  turbulence  the  shock  would  eventually 
. attain  a time-average  asymptotic  thickness.  An  estimate  of  this  asymptotic  thickness  was  made  for  a particular 

case  and  it  gave  a thickening  only  twice  as  great  as  that  produced  by  a single  wave.  The  amplitude  and 
wavelength  were  again  limited  by  computational  restraints  but  it  is  argued  that  the  effects  of  increasing  either 
or  both  would  be  offset  by  the  greatly  increased  tendency  of  the  shock  to  steepen  as  each  successive  wave 
passed  through  it.  It  seems  unlikely,  therefore,  that  the  process  considered  here  could  produce  a thickening  of 
the  magnitude  experienced  by  the  shock  wave  associated  with  the  sonic  boom  on  passing  through  the 
atmosphere. 
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The  problem  is  simulated  in  one-dimension  by  a variation  of  the  The  problem  is  simulated  in  one-dimension  by  a variation  of  the 

piston  problem'  in  which  pressure  waves  representing  turbulence  are  piston  problem'  in  which  pressure  waves  representing  turbulence  are 

made  to  pass  through  a shock  wave.  The  governing  partial  differential  made  to  pass  through  a shock  wave.  The  governing  partial  differential 

equations  are  solved  using  a finite  difference  technique.  equations  are  solved  using  a finite  difference  technique. 

The  results  of  the  calculations  show  that  a shock  wave  tends  to  become  The  results  of  the  calculations  show  that  a shock  wave  tends  to  become 


